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A traditional maritime radar system is utilized for ship detection and 
tracking through onshore transmitters and receivers. However, it faces 
challenges when it comes to detecting small boats. In contrast, unmanned 
surface vehicles (USVs) have been designed to monitor maritime traffic. 
They excel in detecting vessels of various sizes and enhance the capabilities 
and resolution of maritime radar systems. Nevertheless, just like conventional 
radar systems, USVs encounter difficulties due to environmental interference 
and clutter, affecting the accuracy of target signal detection. This research 
proposes a comprehensive numerical assessment to tackle the clutter issue 
associated with USVs. This involves gathering clutter signal data, performing 
numerical analysis, and employing distribution fitting techniques that leverage 
mathematical distributions to unravel data complexity. The root mean square 
error (RMSE) is applied in this analysis to validate the efficacy of the 
distribution model. The results of this study aim to formulate a clutter model 
that can enhance radar performance in detecting small vessels within 
cluttered environments. 
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1. INTRODUCTION 

Conventional radar systems, which are responsible for transmitting and receiving signals, are 
typically placed on land to detect marine vessels [1]. The operation of this system consists of sending 
electromagnetic signals in the direction of the target and detecting and capturing those signals that are 
reflected back to the receiver. Not only does the signal that was received reveal the location of the target, but 
it also makes it possible to conduct an analysis of the target’s size, direction, and velocity [2]. The 
conventional radar setup, on the other hand, frequently has a restricted number of degrees of freedom, which 
can be problematic when attempting to extract information from the station for use in remote sensing 
applications [3]. 

Maritime radar systems employ various frequencies, such as S-band, X-band, K-band, and C-band, 
to cater different requirements [4]. Extensive research substantiates that higher frequencies offer a more 
favorable prospect for detection [5]. By utilizing a 24 GHz frequency, superior detection resolution can be 
achieved, effectively mitigating false alarms. To circumvent delays associated with pulse waves, frequency 
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modulated continuous wave (FMCW) signals can be harnessed within radar systems. Numerous studies have 
successfully implemented a 24 GHz FMCW frequency in their investigations [6]. The utilization of K-band 
frequencies aids in minimizing detection errors [7]. Furthermore, the maritime radar community has 
employed a 24 GHz FMCW frequency for target detection and localization [8]. 

Figure 1 illustrates the propagation of an electromagnetic signal from the shoreline to a sea target 
(black dashed arrow), with subsequent reflection back to the receiver (red dashed arrow). The fidelity of the 
returned signal is compromised by unwanted reflections and refractions caused by factors like sea waves, 
atmospheric conditions, and avian disturbances known as “sea clutter” [9]. Previous works [3], [10], [11], 
focused on radar detection of large vessels such as ferries and cargo ships but challenges persist for smaller 
vessels like fisherman’s boats. Researchers are now exploring unmanned surface vehicles (USVs) as a novel 
approach to improve resolution and capabilities of maritime radar systems [12]. 
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Figure 1. Example of conventional radar 


USVs, autonomous vehicles operating across terrains [13], are proven effective in surveillance and 
monitoring [14]-[16], widely used in civilian and military domains. In radar detection, bistatic radar 
implementation, eliminating bulky machinery [3], has promising results. The proposed approach involves 
leveraging USVs by equipping them with radar sensors and strategically deploying them at sea locations. 
However, sea clutter still affects detection probability, especially for small vessels. 

Danklmayer et al. [17] introduces a novel sea-to-sea wave propagation method to mitigate sea 
clutter in maritime radar applications. This approach enhances sensing signals, improving radar cross section 
(RCS) returns for small boats through space diversity effects [18]. The study validates a newly developed 
clutter model through comprehensive numerical evaluation, aiming to improve existing radar systems and 
reduce clutter. However, literature on sea-to-sea wave propagation lacks numerical evaluation models for a 
comprehensive description of radar clutter in such systems. 

Signal transmission from shore to sea creates triangulated reflections due to sea waves, resulting in 
disruptive sea clutter spikes [19]. To address this challenge, the study introduces a novel detection system 
utilizing a sensor mounted on a USV. The system gathers and analyzes data to investigate sea clutter 
characteristics, with a focus on sea-to-sea wave propagation and the utilization of space diversity. 

The study aims to address adverse sea clutter effects on signal quality and enhance overall detection 
performance. This involves determining a suitable distribution that accurately represents clutter characteristics, 
utilizing a 24 GHz radar system for improved resolution and discrimination capabilities. The paper begins with 
a detailed explanation of radar requirements derived from a Fresnel zone study, followed by an examination of 
experimental configuration and numerical analysis of clutter phenomena. The study includes a comprehensive 
analysis of the distribution model. The final section outlines future work for advancing a maritime clutter model 
tailored for US Vs. Experimental results can improve radar system performance for small vessel detection and 
weather forecasting through numerical evaluation and data fitting [20]-[22]. 


2. METHODOLOGY 

The comprehensive study as depicted in Figure 2, begins with the determination of the important radar 
requirements for the placement and configuration of the transceiver. This is achieved by utilizing the concept of 
the Fresnel zone. Following this, the experimental configuration is customized to align with the capabilities of 
the Pusat Penyelidikan Air Malaysia (NAHRIM) facility, facilitating the production of artificial sea waves with 
predetermined heights and intervals. After a collection of cluttered data, the Hilbert transform is utilized to 
extract the amplitude, thereby improving the quality of the data. The processed data is utilized to reveal the 
hidden Doppler information signal. By employing the parameters obtained from the analysis above, a thorough 
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evaluation of clutter measurement is conducted by constructing probability distribution functions (PDFs) that 
encompass log-normal, log-logistic, gamma, and weibull distributions. In this study, various quantitative 
metrics are employed to evaluate the distributions under investigation. These metrics include maximum 
likelihood estimation (MLE), standard deviation (STD), goodness-of-fit (GOF), and root mean square error 
(RMSE). The computed values of these metrics for each distribution are then presented and discussed in the 


subsequent results section. 


RADAR MEASUREMENT DATA SIGNAL DATA 
REQUIREMENT SETUP COLLECTION PROCESSING COLLECTION 
+ Fresnel zone e Water level e Raw data e Hilbert transform * PDF 
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configuration e Devices * GOF 
e RMSE 


Figure 2. Clutter evaluation process flow 


2.1. Radar requirement estimation 

Mounting the transceiver in radar requirement estimation involves considering the Fresnel zone for 
optimal signal propagation, crucial for detecting surrounding objects and minimizing blind spots. A 
comprehensive monitoring system allows a 360-degree view of the radar’s surroundings [23]. The height of 
the radar system’s antenna directly affects transmitting and receiving antenna angles, as extensively studied 
at a 5 GHz frequency [4]. These investigations offer valuable insights into radar system design, ensuring 
optimal functionality and performance [24]. 

Depending on the size of the USV, placing transmitting and receiving antennas at a height of 1 m 
above the ground is sufficient to eliminate destructive ground waves that could interfere with the radar 
processing. Figure 3 is the illustration of the concentric ellipse of transmission of the antenna. The total 
distance from the transmitting and receiving antenna is 10 m which d; and dz are both 5 m while h, and h, are 
both 1 m. The Fresnel zone can be calculated by using the basic (2): 


Radius (mts.) = 17.31 x Es o 
4 xf (in GHz) 


where D is the distance of the transmitting antenna to the center of the transmission path, and r is the radius 
of the concentric ellipse respectively as shown in Figure 3. 


2— |i 
RQ AAAS SS 
D 


Figure 3. Concentric ellipse of transmission 


Utilizing a distance, D of 10 m and an operating frequency of 24 GHz, the radius, r of the Fresnel 
zone can be determined where the calculated radius is found to be 0.177 m. Notably, this value represents the 
minimum height requirement for the antenna to prevent ground reflections and maintain optimal 
performance. Failure to adhere to this minimum antenna height would result in undesired beamwidth 
bouncing off the ground, leading to undesirable ground reflections and subsequent impact on the radar 
system's overall performance. Therefore, ensuring a minimum antenna height of 0.177 m is essential to 
mitigate these adverse effects and maintain the desired operational characteristics of the radar system. This 
calculation and requirement contribute to the accurate design and implementation of radar systems, 
optimizing their performance in practical applications. The height of antenna placement in this study is 
placed 1 and 1.5 m above sea water level which exceeds the minimum height of antenna placement as per 
calculated above. 
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2.2. Measurement setup 

Data collection and clutter measurement for this study take place at NAHRIM, a renowned research 
institution with expertise in water-related investigations. To ensure accuracy, experimental parameters are 
carefully set and maintained in a controlled environment. Due to limitations in wave height at NAHRIM, a 
1:100 scale-down approach is employed, allowing for accurate representation of real-world conditions while 
facilitating precise measurements and data collection. 

The experiment is conducted in a basin of 30 m?. An artificial sea wave is generated by using a wave 
generator at the end of the basin and the transceiver is mounted facing toward the wave generator as shown in 
Figure 4(a). During the measurement, an OmniPreSense radar sensor is used to collect the clutter signals where 
the wave signal is propagated on the water surface, bounced back, and received by the receiver antenna. Each 
data is collected for 60 s and the total of data collected for this analysis is 30 data. The collected signal includes 
noise and clutter along the parameter. It is initially processed to determine the average values, which are then 
subtracted from the raw data to obtain the mean value. Subsequently, the processed data underwent a Hilbert 
transformation to extract the amplitude of the clutter. The resulting processed data, characterized by its 
enhanced signal quality, will be utilized to determine the enveloped signal of the Doppler information. 

Figure 4(b) presents a side view of an illustrative depiction of a transceiver system positioned above the 
water surface, highlighting the key elements and processes involved in sea-to-sea wave propagation. The 
transceiver antenna, denoted as TX (transmitter) and RX (receiver), serves as the source and destination for signal 
transmission, respectively. Upon transmission from the transceiver antenna, the electromagnetic signal propagates 
across the water surface, exhibiting a characteristic bouncing effect, before ultimately returning to the transceiver 
antenna. This phenomenon, known as sea-to-sea wave propagation, plays a pivotal role in the formation and 
behavior of the transmitted signal. Notably, the height of the transceiver antenna is represented as h, while Iw and 
hw symbolize the wave interval and wave height, respectively. The blue and red arrows in Figure 4 visually 
represent the propagation of the signal along the water’s surface. A comprehensive understanding of these 
parameters and the underlying sea-to-sea wave propagation mechanisms is crucial for accurately assessing the 
impact of sea clutter on signal quality and devising effective mitigation strategies. The transceiver configuration 
during data collection is shown in Figure 5 where the antenna is securely mounted on a pole positioned above the 
seawater surface before the wave is generated, while the parameters used are tabulated in Table 1. 


Tx /Rx 


(a) (b) 


Table 1. Measurement 
parameters and value 


Parameters Value 

Antenna land1.5m 

height (h) 

Antenna Transceiver 

configuration 

Wave height 10 and 20 

(hw) cm 

Wave 2s 
Figure 5. Actual transceiver configuration interval (Iw) 


2.3. Numerical analysis of clutter data 

This section delves into the mathematical models utilized to represent distinct clutter signal samples. 
The determination of parameters for each model is based on mathematical equations specifically tailored for this 
purpose. The core objective of this analysis is to ascertain the parameter values of the optimal-fit distribution 
model by leveraging the measured clutter data. The primary focus of this analysis centers around the PDF of 
Weibull, Gamma, log-logistic, and log-normal distributions as summarized in Table 2, which have garnered 
considerable interest despite the presence of several alternative models in the existing literature. 
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Table 2. Distribution model 
Description 


No. Distributions Mathematical format 


1 Log-normal The log-normal distribution, defined by u and o, models positive data by (f(x|u, 0) = 
distribution exponentiating a normally distributed random variable, with u as the 1 exp { -0nx- oh; 
location and o as the scale. x27 20° i 
x 20 
2 Log-logistic The log-logistic distribution, similar to the log-normal distribution, Falu od = at ee ; 
distribution features heavier tails. It is often employed in survival analysis for ox (1467)? 
parametric modeling of events with increasing and decreasing rates. i z0 
where z = 28- 
O 
3 Gamma The Gamma distribution model is a continuous distribution where it is the (f(xla, b) = 1 ya-1 ar 
distribution sum of exponentially distributed random variables. b*Ia) 
4 Weibull Is a continuous probability distribution using positive only values of data. (f(xla, b) = b Ör: eG” 
distribution Strictly positive values of the shape parameter b and scale parameter a. aa 


The numerical evaluation of the clutter data will go through all these distributions and calculations 
such as MLE, STD, GOF, and RMSE are made for each distribution model. To identify and validate the 
model that exhibits the least error between the clutter data and the statistical model, a GOF test approach 
employing the RMSE method is employed in this analysis. 

MLE is utilized to determine the parameters of fitting for each distribution. The parameters include 
the values of (u) and (o), which are used with the log-normal and log-logistic distributions. In contrast, alpha 
(a) and ($) values are determined for the Gamma and Weibull distributions. STD is the square root of the 
variance data sample to determine an error for each distribution parameter, as shown in (2), where p for 
population data and ois sample data. 


Sih (2) 

The GOF for the statistical models is evaluated using the RMSE method [25]. RMSE, defined as the 
standard deviation of the residuals or prediction errors, serves as a measure of the accuracy of the fitted 
distributions. The distribution that exhibits the lowest percentage error is considered the most suitable fit. To 
assess the GOF, RMSE is calculated for each distribution by comparing the amplitude values of the sample 
clutter data (i) with the corresponding amplitude values of the statistical model (%), The RMSE value is 
computed using (3), where i represents the sample clutter data number. This rigorous evaluation contributes 
to the assessment of the GOF and enhances the reliability of the statistical modeling process. 


RMSE= [+ n (G 4) (3) 


3. RESULTS AND DISCUSSION 

The findings of this analysis have been divided into two sections: clutter processing and numerical 
analysis of a clutter signal that represents a distinct wave height corresponding to a specific antenna height. 
For numerical analysis, the following distribution models are chosen for each category of wave height: 
Weibull, Gamma, log-logistic, and log-normal. Table 3 contains tabular representations of the calculated 
parameters, STD, and RMSE values for all models, respectively. 


oO = 


Table 3. Distribution model parameters, STD, and percentage of RMSE 


Parameters 1D And S EDE : 
STD Err Log-logistic Log-normal Gamma Weibull 
Im 10 cm STD u=2.0723 u=1.9872 a=2.9006 a=9.8326 
antenna wave 6.8810 o=0.3574 o=0.7226 B=3.0176 B=2.0157 
height STD Err RMSE, =0.0034 RMSE,;n=0.0040 RMSEg=0.0033 RMSEw=0.0030 
0.0089 %err.=0.3426 %err .=0.4024 “err .=0.3318 Y%err .=0.3010 
20 cm STD u=2.1191 u=2.0377 a=3.1806 a=10.1741 
wave 6.4416 o=0.3336 o=0.6838 B=2.8461 B=2.0916 
height STD Err RMSE, =0.0044 RMSE;n=0.0052 RMSEg=0.0044 RMSEw=0.0040 
0.0083 %err .=0.4352 %err .=0.5160 %err .=0.4364 Yerr .=0.3997 
1.5m 10cm STD u=2.4030 u=2.3507 a=4.1578 a=13.4062 
antenna wave 8.3440 o=0.2942 o=0.5513 B=2.8598 B=2.3662 
height STD Err RMSE, =0.0034 RMSE,n=0.0037 RMSEg=0.0033 RMSEw=0.003 1 
0.0108 “err =0.3393 %err =0.3682 %err =0.3277 %err =0.3108 
20 cm STD u=2.1583 u=2.0813 a=2.6761 a=10.9842 
wave 8.1620 o=0.3813 o=0.7285 B=3.6520 B=1.8715 
height STD Err RMSE, =0.0033 RMSE;n=0.0036 RMSEg=0.0030 RMSEw=0.0027 
0.0106 Yerr =0.3294 %err =0.3637 %err .=0.2958 %err =0.2681 
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3.1. Clutter processing 

As mentioned earlier, a comparative analysis is conducted in the time domain to evaluate the 
processed data obtained from antennas positioned at heights of 1 and 1.5 m. Throughout the experiments, the 
wave height is kept constant (either 10 cm or 20 cm) while the antenna heights are varied. The processed 
time domain clutter signal, focusing on the initial 60 s, is displayed in Figures 6(a) and (b) for antenna 
heights of 1 m and 1.5 m, respectively, while considering a wave height of 10 cm. In Figure 6(a), 
corresponding to the | m antenna height, the amplitude clutter signal is observed to range from 12 to 30, with 
an average amplitude value of 18. Conversely, Figure 6(b) portrays the clutter signal for the 1.5 m antenna 
height, exhibiting an amplitude range of 22 to 38, and an average amplitude value of 27. 


, Amplitude(Volt) 

; aak 
Amplitude (Volt) 

; Amplitude(Volt) 


0 10 20 30 40 50 60 o 10 2 30 40 50 60 0 10 20 30 40 50 60 0 10 20 30 40 50 60 
Time (second) Time (second) Time (second) Time (second) 


(a) (b) 


Figure 6. Time domain clutter signal for with 1 m and 1.5 m antenna height; (a) 10 cm wave height and 
(b) 20 cm wave height 


Further investigation is conducted to examine the impact of a higher wave height, specifically set at 
20 cm. Figure 6 illustrates the time domain clutter signal obtained under two different antenna heights, 
namely | m (Figure 6(a)) and 1.5 m (Figure 6(b)). In the case of an antenna height of | m, the range of the 
amplitude clutter signal, as depicted in Figure 6(a), varied from 17 to 40. Notably, the average value of the 
amplitude signal is measured to be 22. On the other hand, when the antenna height is increased to 1.5 m, the 
range of the amplitude clutter signal exhibited a slightly different pattern, ranging from 15 to 40. 
Remarkably, the average value of the amplitude signal in this scenario is recorded as 24, as illustrated in 
Figure 6(b). This analysis provides valuable insights into the influence of wave height and antenna height on 
the characteristics of the clutter signal, contributing to a deeper understanding of radar system performance. 

Figure 7(a) depicts the amplitude-enveloped signal for 20 cm wave heights, which exhibits 
significantly larger magnitudes than the signal for 10 cm wave heights. This variation in amplitude can be 
attributed to the spatial limitations imposed by lower antenna heights, which result in reduced surface 
coverage. Despite these limitations, it is notable that the radar system collects a greater amount of clutter in 
the case of higher wave heights (20 cm versus 10 cm). This finding highlights the impact of varying sea wave 
heights on system noise, indicating that higher waves generate more significant interference in radar clutter. 
Meanwhile, Figure 7(b) depicts a comparison of the enveloped signals of 10 cm and 20 cm wave heights for 
a fixed 1.5 m antenna, revealing that the amplitudes for each wave height are only slightly different. The 
amplitude of interference is not overly significant since larger antennas cover a larger area and therefore 
capture more interference into the system. 


65 65 p, 
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60 ---20 em| 60 ---20 cm 
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Figure 7. Envelope signal; (a) 1 m antenna height and (b) 1.5 m antenna height 
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3.2. Distribution fitting and numerical analysis of clutter data 

Figure 8 presents a comprehensive analysis of the clutter signal probability plot, showcasing a 
comparison of four fitting models: log-normal, log-logistic, Gamma, and Weibull distributions. The 
probability plots are displayed for both a wave height of 10 cm and 20 cm, considering antenna heights of 
1 mand 1.5 m. In Figure 8(a), the probability plot for the clutter signal with a 10 cm wave height is depicted 
for antenna heights of 1 and 1.5 m, respectively. Additionally, Figure 8(b) represents the probability plot for 
the clutter signal with a wave height of 20 cm, again for antenna heights of 1 and 1.5 m. 


0.05 - 0.05 0.05 0.05; — 
Clutter Clutter Clutter | Clutter 
Log-Normal Log-Normal Log-Normal | Log-Normal 
0.04 + Log-Logistic 0.04 Log-Logistic 0.04 + Log-Logistic | 0.04 Log-Logistic 
+ Gamma * Gamma < Gamma | + Gamma 
2003 | > Weibull 20.03 Weibull 2003 Weibull | 2003 Weibull 
2 a a a 
f f f 
à 0.02 È £ 20.02 
ra 
0.01 0.014 > \ 
0 ; ; À of 
0 10 20 30 40 50 0 10 20 30 40 50 0 10 20 30 40 50 0 10 20 30 40 50 
Clutter Clutter Clutter Clutter 
(a) (b) 


Figure 8. Distribution fitting for with 1 m and 1.5 m antenna height; (a) 10 cm wave height and (b) 20 cm 
wave height 


Remarkably, the Weibull distribution emerged as the most suitable fitting model across all 
probability plots. This implies that the Weibull distribution effectively captures the statistical characteristics 
of the clutter signal. Furthermore, it is noteworthy that the probability plots exhibit greater scatter for the 
higher wave height scenario. This observation can be attributed to the increased clutter signal generated by 
higher wave heights, which consequently exerts a significant impact on radar system performance. These 
findings provide valuable insights into the selection of appropriate distribution models and the understanding 
of radar performance in the presence of varying wave heights. 


4. CONCLUSION 

Radar-based USVs can overcome the limitations of conventional maritime radar in detecting small 
vessels. However, it still suffers from the effect of clutter, which can reduce the probability of detection; this 
can be resolved by analyzing the optimal distribution and comprehending the characteristics of clutter. The 
Weibull distribution is the best-fitting distribution in this investigation, based on the preceding finding. This 
is confirmed by the RMSE value, which provides the lowest error rate, and this evaluation promises a higher 
resolution for vessel detection in maritime radar. These results also demonstrate a significant differential 
effect, suggesting that antenna placement at a higher altitude above sea level results in increased surface 
coverage. However, it captures a greater amount of noise and interference from the environment. As a result, 
it is crucial to take antenna height into account when designing and deploying radar systems, as it has a direct 
impact on the characteristics of clutter and the overall efficacy of the radar system. This finding will 
contribute to the modeling of maritime USV clutter, which can be used to reduce false alarms and enhance 
the detection resolution of small vessels. 
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